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ABSTRACT 


A new method is described for determining which stars in a navigation star 
catalog are near a star sensor field-of-view (FOV). This method assumes that 
an estimate of spacecraft inertial attitude is known. Vector component ranges 
for the star sensor FOV are computed, so that stars whose vector components lie 
within these ranges are near the star sensor FOV. This method requires no 
presorting of the navigation star catalog, and is more efficient than 
traditional methods. 
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1.0 INTRODUCTION 


Many spacecraft maintain inertial attitude knowledge by sighting navigation 
stars with star sensors. The process of determining exactly which navigation 
stars are in a star sensor field-of-view (FOV), given an estimate of spacecraft 
attitude, is computer intensive. Consequently, spacecraft first determine 
which navigation stars are near the star sensor FOV, and then determine 
which nearby stars are in the star sensor FOV. This paper presents a new 
method of determining which navigation stars are near a star sensor FOV that 
is more efficient than previous methods. 


2.0 BACKGROUND 

There are two classes of methods for determining which navigation stars are 
near a star sensor FOV: those methods that do not require presorting the 

navigation star catalog and those methods that do require presorting the 
navigation star catalog. 

The primary no-presorting method is the scalar-product method. This method 
is used on the Space Shuttle. It takes advantage of the simple scalar-product 
relation 


0 = cos'^b • s) 

where b is a unit star sensor boresight vector in inertial coordinates, s is a 
unit navigation star position vector in inertial coordinates, and 0 is the angle 
between b and s. If 0 is less than the angular radius of the star sensor FOV, 
then s is near the star sensor FOV. Obviously, the scalar-product method is 
simple, but it is also computer intensive. 

There are numerous presorting methods that determine which navigation 
stars are near a star sensor FOV. Although these methods are normally less 
computer intensive than the scalar-product method, they do introduce some 
additional problems. For example, some algorithms require that the 
navigation star catalog be presorted into predefined zones on the celestial 
sphere. These algorithms are complicated and require additional computer 
data storage to index the navigation stars. 

It is desirable to have a method for determining which navigation stars are 
near a star sensor FOV that is simple, fast, and requires no presorting of the 
navigation star catalog. The method described below satisfies these 
requirements: it is a no-presorting method that is faster than the scalar- 

product method. 
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3.0 METHOD 


Let us assume that we have a star sensor with a known FOV and an associated 
navigation star catalog. In addition, let x, y, and z be the principal axes of the 
inertial coordinate system in which the coordinates of the navigation stars are 
cataloged. Then we may define orthogonal bands on the celestial sphere that 
are normal to the x, y, and z axes; and whose unique intersection I bounds a 
circular cap on the celestial sphere that is circumscribed about the star sensor 
FOV. This concept is illustrated in Figure 1 below. 
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Figure 1 - Intersection of orthogonal bands containing star sensor FOV. 

To determine which navigation stars are near the star sensor FOV, we need 
simply determine which navigation stars are in £. Since the orthogonal bands 
illustrated above are normal to the x, y, and z axes, £ may be defined by its 
projections onto x, y, and z. Let the intervals [x m in> x m axl. [ymin* ymaxl- an ^ 
[zmin. Zmaxl be the projections of £ onto x, y, and z respectively. Then to 
determine which navigation stars are in £, and thus near the star sensor FOV, 
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we need simply search the navigation star catalog for stars whose inertial 
coordinates are contained in [x m j n , x max ], [ymin> Ymaxl' an ^ t z min> z maxl- 


4.0 ALGORITHM 

To compute the coordinate interval [x m i n , x max ], let b = (b x , by, b z ) be the star 
sensor boresight vector in inertial coordinates, a x be the angle between the 
boresight vector and the x-axis, and P be the diagonal radius of the star sensor 
FOV. These parameters are illustrated in Figure 2 below. 
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Figure 2 - Star sensor FOV within coordinate interval [x m i n , x max ] . 
The angle between the boresight vector and the x-axis, a x , is given by 

a x = cos’^bx), where 0 < a x < 7t. 
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The coordinate interval [x m i n , x max ] may be computed as follows: 

if (a x - P) < 0 then 
x max = 1 

Xmin = cos(a x + P) 

elseif (a x + P) > n then 
Xmax = cos(a x -P) 

x min = " 1 

else 

xmax = cos(a x - P) 

xmin = cos(a x + P) 

endif 

The intervals [y m in> ymaxl and [z m i n , z m axl can be computed in a completely 
analogous manner. 

Once the coordinate intervals [x m i n . x m axl. [ymin. ymaxl. and [z m i n , z max ] have 
been computed, we must search through the navigation star catalog for stars 
whose inertial coordinates are contained in these intervals, and thus in L and 
near the star sensor FOV. The following is an example of how this process may 
be implemented. 

loop Star Catalog 

if s x > x m i n then 
if s x < Xmax then 
if Sy > ymin then 
if Sy < ymax then 

if Sz > z m j n then 

if Sz < z max then 

Star is near star sensor FOV 
endif 
endif 
endif 
endif 
endif 
endif 
end loop 

The algorithm above may be improved by ordering the "if' statements so that 
the least-probable "if' statements appear first. 
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5.0 TEST RESULTS 


The results of a comparison between the scalar-product method and the new 
method described above are given in Table 1 and Table 2 below. The tests were 
run on a Macintosh IIx using Absoft MacFortran/020, Version 2.4. A constant 
star sensor FOV diagonal radius of 8 degrees was assumed. Table 1 summarizes 
data generated using a math coprocessor, and Table 2 summarizes data 
generated without using a math coprocessor. 

Table 1. Comparison of the scalar-product method and the new method using a 
math coprocessor. 


Number 
of Trials 

Number 
of Stars 


New Method 
(sec) 

Ratio of Times 
(S-P : New) 

1000 

100 

4.72 

1.93 

2.4 : 1 

1000 

1000 

46.85 

16.33 

2.9 : 1 

1000 

5000 

234.18 

80.05 

2.9 : 1 


Table 2. Comparison of the scalar-product method and the new method without 
using a math coprocessor. 


Number 
of Trials 

Number 
of Stars 

Scalar-Product 
Method (sec) 

New Method 
(sec) 

Ratio of Times 
(S-P : New) 

1000 

100 

27.33 

14.37 

1.9 : 1 

1000 

1000 

273.17 

26.17 

10.4 : 1 

1000 

5000 

1365.85 

78.73 

17.3 : 1 


The test results indicate the new method is superior to the scalar-product 
method. Improvement was shown over the scalar-product method, both with 
and without a math coprocessor. 
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6.0 CONCLUSIONS 


A new method for quickly identifying which stars in a navigation star catalog 
are near a star sensor FOV has just been described. This method is faster than 
the traditional scalar-product method, yet retains the advantage of requiring 
no presorting of the navigation star catalog. 

If a faster algorithm is desired, the navigation star catalog may be presorted 
by unit vector components. This will increase computational speed, but it is 
more complicated and will require additional computer data storage to index 
the stars. 

The method is presented as an application for onboard spacecraft attitude 
determination, but it may also apply to other areas such as astronomy. 
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